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EDITOR'S  PREFACE 

One  of  the  more  frequently  quoted  papers  in  fluid  dynamics  is  the 
article  by  Karl  Pohlhausen,  “Zur  naherungsweisen  Integration  der 
Differentialgleichungen  der  laminaren  Grenzschicht,”  which  appeared 
on  pages  252-268  of  Volume  I  of  the  Zeitschrift  fur  angewandte 
Mathematik  und  Mechanik.  Consequently,  it  seems  appropriate  to 
make  an  English  translation  of  the  original  German  paper  widely 
available. 

While  the  approximate  integration  technique  which  the  article 
outlines  has  been  superseded  in  large  measure  by  substantial  im¬ 
provements,  the  kernel  of  all  useful  methods  is  given  in  the  original 
paper;  in  fact,  while  preparing  a  catalogue  of  integral  equations,  I 
was  astounded  to  record  the  number  of  variations  on  Pohlhausen’s 
method  that  have  been  proposed.  Moreover,  many  other  important 
results  of  permanent  value  to  fluid  dynamics  are  presented,  e.g.,  a 
mathematical  derivation  of  K<irm>in’s  momentum  equation  from  the 
boundary  layer  equation  and  the  mathematical  description  in  finite 
closed  form  of  the  steady  two  dimensional  laminar  How  in  a  con¬ 
verging  channel  according  to  boundary  layer  theory. 

So  much  for  the  raison  d’etre,  the  translation  is  the  work  of 
Roland  C.  Anderson  and  represents  a  scholarly  effort  that  was  sub¬ 
mitted  to  the  Graduate  School  of  the  University  of  irida  in  lieu 
of  an  examination  on  a  reading  knowledge  of  German  it  is  one  re¬ 
quirement  for  a  doctorate  in  aerospace  engineering.  Richard  L. 
Fearn  has  checked  the  entire  paper  in  minute  detail  and  has  made 
numerous  minor  corrections.  I  hope  that  he  caught  all  of  the  errors; 
if  not,  please  let  me  know.  The  drawings  are  the  expert  works  of 
Karl  P.  Rhodes  and  Roland  C.  Anderson. 

Finally,  an  appropriate  acknowledgment  for  the  gracious  financial 
and  psychological  support  of  the  Office  of  Naval  Research  under  Con¬ 
tract  Nonr-580(14)  is  a  debt  that  I  now  discharge  with  sincere  g  ati- 
tude  and  great  pleasure. 

Knox  Millsaps 

Gainesville,  Florida 
15  August  65 
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IV 


THE  APPROXIMATE  INTEGRATION  OF  THE 
DIFFERENTIAL  EQUATION  FOR  THE 
LAMINAR  BOUNDARY  LAYER 


by 

Karl  Pohlhausen 

The  extraordinary  mathematical  difficulties  which  occur  in  the 
integration  of  the  differential  equations  for  fluid  motion,  especially 
with  the  added  complications  due  to  viscosity,  and  which  are  due  to 
the  essential  nonlinear  character  of  the  equations  has  caused  a  dicho¬ 
tomy  between  “Hydrodynamics”  and  “Hydraulics”  in  the  develop¬ 
ment  of  the  theory  of  fluid  flow.1  On  one  hand,  researchers  neglected 
viscous  effects  in  order  to  obtain  simple  equations  and  achieved  rigor 
with  substantial  discrepancies  between  the  calculated  and  observed 
Hows.  On  the  other  hand,  a  technology  which  was  necessitated  by  a 
consideration  of  the  actual  behavior  of  a  fluid  developed  into  hydrau¬ 
lics  which  was  a  separate  science  of  fluid  motion.  In  hydraulics,  the 
exact  equations  were  replaced  by  empirical  assumptions  and  intuitive 
considerations  whose  results  agree  essentially  with  reality. 

In  recent  times,  researchers  have  tried  to  bring  hydrodynamics 
and  hydraulics  into  agreement  by  including  viscosity  in  the  mathe¬ 
matical  theory  and  by  looking  deeper  into  hydraulics.  The  first 
general  beginnings  which  included  frictional  considerations  were 
made  by  Stokes.2  He  was  able  to  calculate  the  resistance  of  a  sphere 
moving  in  a  fluid  with  very  larg.  viscosity.  The  assumption  behind 
this  solution  was  that  the  motion  will  essentially  be  determined  by 
viscous  effects,  thus  either  the  viscosity  is  very  large  or  the  product 
of  the  body  dimension  and  the  velocity  is  small.  A  systematic 
simplification  tor  the  flow  of  a  fluid  of  small  viscosity,  given  first  by 


iSee  l*h.  Forchheimer,  Hydraulik,  Leipzig.  1914.  p.  70. 

jG.  Stokes,  On  the  effect  of  the  internal  friction  of  fluids  on  the  motion  of 
pendulums,  Prof.  Camb.  Ml.  Sac.,  9,  1851  \. Men .  Paprrs.  },  p.  I). 
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Prandtl,3  led  to  the  theory  of  the  “boundary  layer"  which  was  ex¬ 
amined  in  detail  and  applied  to  selected  examples  by  his  students 
Blasius,4  Boltze5  and  Hiemenz.8 

An  approximate  method  of  integration  for  the  differential  equation 
of  the  boundary  layer  will  be  developed  in  the  following,  and  its 
usefulness  will  be  shown  by  examples. 

]. Derivation  of  the  Differential  Equation  for  the  Boundary  Layer. 

Prandtl  considered  the  motion  of  a  fluid  about  a  fixed  body  sub¬ 
mersed  in  a  stream  or,  equivalently,  about  a  body  moving  at  a  con¬ 
stant  velocity  in  a  still  fluid.  He  arrived  at  the  differential  equation 
for  the  boundary  layer  by  assuming  the  coefficient  of  viscosity  to  be 
very  small.  At  a  large  distance  from  the  body  the  friction  will  have 
no  influence;  consequently,  the  motion  far  from  the  body  will  be 
potential  flow.  Only  in  a  very  thin  layer  near  the  body,  “the  boundary 
layer,”  will  there  be  any  deviation  from  potential  flow,  and  this 
deviation  is  caused  by  the  fluid  clinging  to  the  wall.  If  one  denotes 
the  thickness  of  the  boundary  layer  by  8,  the  order  of  magnitude  of 
each  term  in  the  exact  differential  equations  may  be  estimated.  As 
Blasius  has  rigorously  shown,  the  neglect  of  terms  of  small  order 
of  magnitude  led  to  the  Prandtl  boundary  layer  equations. 

We  intend  to  obtain  the  same  differential  equation  here  in  an¬ 
other  way,  i.e.,  by  asymptotic  expansion.7 

We  start  with  the  Navier-Stokes  equations  which  describe  the 
flow  of  a  fluid  with  friction,  and  limit  ourselves  to  steady  flow  in  a 
plane.  First  of  all,  we  shall  carry  out  the  calculation  for  the  flow 


3L.  Prandtl,  Verhand.  des  III.  Intern.  Math.  Kongres,  1904  (Heidelberg),  Leip¬ 
zig  1905,  p.  484.  An  English  translation  is  available  as  NACA  Technical  Memoran¬ 
dum.  No.  452,  1928. 

*H.  Blasius,  Dissertation  GOttingen,  1907.  Zeitschr.  filr  Math.  u.  Phys„  56,  1908. 
p.  1.  An  English  translation  is  available  as  NACA  Technical  Memorandum,  No. 
1256,  1950. 

»E.  Boltze,  Dissertation  Gdttingen,  1908. 

<'K.  Hiemenz,  Dissertation  GOttingen,  1911.  Dinglers  Polyt.  Journal,  526,  1911, 
p.  521. 

’The  derivation  is  universal  in  agreement  with  that  in  the  preceding  work  by 
von  Kdrmdn  which  is  available  in  English  as  NACA  Technical  Memorandum,  No. 
1092,  1946;  on  these  grounds  repetition  is  perhaps  not  inappropriate. 
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along  a  wall  at  y  =  0.  The  equations  are  well  known  to  be 

1  3p 

3x 


3u  3u 

u - f-  V - 

3y 

3v  3v 

u - j-  V -  : 

3x  3y 


p  3x 
1  3p 

p  3y 


•  -(-  vAu 


+  vAv 


3u  3v 

3x  3y 


(1) 


Moreover,  if  the  fluid  must  be  absorbed  on  the  wall  of  the  body, 
it  must  be  true  also  that  u  =  0  and  v  =  0  at  y  =  0. 

The  continuity  equation  may  be  integrated  by  the  introduction  of 
the  stream  function  in  such  a  way  that 

3i/r  dtp 

u  = - and  v  = - 

3y  3x 

and  we  obtain  then  from  the  first  two  equations  in  (1) 
dtp  3  (At p)  dtp  3  (Atp) 


3y 


3x 


3x 


3y 


=  vAAtp 


(2) 


We  introduce  dimensionless  values  by  dividing  all  lengths  by  a 
characteristic  length  of  the  geometry,  a,  and  the  velocities  by  a  char¬ 
acteristic  velocity,  u  (e.g.,  the  undisturbed  velocity  of  the  fluid  at 
infinity),  and  we  write 

tp  =  aQtp'  (x'.y') 

where  x'  =  x/a  and  y'  =  y/a 

From  this,  we  change  (2)  into  the  form 

2  r  dtp'  3  (A tp')  dtp'  d  (A f  )■ 


Of) 


“1  ad 


3y*  3x'  3x' 

or  with  R  =  aQ/v  (Reynolds  number) 

dtp'  d(Atp')  dtp'  3  (A tp') 


3y'  J 


.  pAAtp' 


3y'  3x' 


3x'  3y' 


-  A  Atp' 


(3) 


With  the  restriction  of  large  Reynolds  numbers  we  require  the  solu- 
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tion  to  be  in  the  form 


t/f'  =  f  (x',y')  +  R8F  (x',z) 
z  =  y'R" 

where  0  and  a  are,  for  the  present,  undetermined  exponents,  f  (x‘,  y') 
must  represent  the  potential  solution  which,  of  course,  is  a  rigorous 
solution  of  (2);  thus, 

Af  <x',y')  =  0 

We  now  consider 

32F  32F  32F  32F 

AF  = - 1 - = - U  R2<7 - 

3x'2  3y'2  3x'2  3z2 


Under  the  assumption  that  R  is  very  large  and  a  is  positive, 
32F/3x'2  can  be  neglected.  Likewise  in 

34F  34F  34F 

AAF  = - f-  2R2" - j-  R4<r - 

3x'4  3x'23z2  3z4 


only  the  last  of  the  three  terms  on  the  right  are  retained.  If  these 
expressions  are  inserted  into  (3),  we  get 

3f  33F  3f  3:,F  z 


R0  +  2 


3y'  3z23x' 

3F  33F 


3z  3x'3z2 


_ J^S  +  217 


3x'  3z3  y' 


+ 


3F  33F 


3x'  3z3 


_R@  +  4 


34F 


3z4 


(4) 


With  increasing  R,  the  right  and  left  sides  of  this  equation  must  go 
uniformly  to  infinity;  thus,  the  restrictions  on  ft  and  a  are 

/3  -|-  =  2/3  +  3a  =  (3  +4<r  -  1 

From  this,  it  is  determined  that 

a  =  V*  and  /3  =  —Vt 

If  we  expand  f  (x',y*)  in  a  Taylor  series  in  y', 

f(x\y')  =  f,  +  f.  (x')y'  +  . . . 

where  f0  is  a  constant,  then  the  wall  of  the  body  must  be  a  streamline, 
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and  the  differential  equation  (4)  becomes 


r 

r . 

&l,  - 

i! 

-  aF 

s*r 

3F 

33F  “ 

34F 

3r 

_  3z 

3x‘3z2 

3x* 

3z3  _ 

3z4 

IU* 


/  3F  \  33  (F  +  f,z)  /  3f,  3F  \ 

\  ^  3z  /  3x'3z2  \  3x'  3x'  / 


3‘(F  +  f,z)  34  (F  +  f,z) 


3z3 


3z4 


as  one  can  easily  prove  by  differentiation.  Letting  (F  -|-  f,z)  equal 
G,  we  obtain 


3G  3:,G  3G  3*G  34G 

3z  3x'3z2  3x'  3z3  3z4 


(5) 


Thus,  in  terms  of  the  function  G,  the  solution  now  reads 
V  =  f  (x'.y')  +  (G  -  f,z)R-v* 
where  the  boundary  conditions  for  G  are 

1)  for  z  =  0  :  G  =  0 

2)  for  z  =  0  :  3G/3z  ==  0 

3)  for  z  =  «  :  3G/3z  =  f, 


We  have,  in  fact,  obtained  a  second  approximation  for  the  flow 
of  a  fluid  under  the  assumption  of  large  Reynolds  numbers,  since  in 
the  vicinity  of  the  wall  if/'  =  GR-,/2,  then  on  the  wall  f(x',  y')  be¬ 
haves  like  f,y'  or  f,z  R“1/2,  and  the  no  slip  condition  is  satisfied  at 
y  =  0.  At  infinity,  G  behaves  like  f,z;  therefore,  G  —  f,z  =  0  leaving 
only  \fi'  =  f  (x'.y'),  i.e.,  the  potential  solution.  If  we  integrate  (5) 
with  respect  to  z,  then 

3G  32G  3G  3-G  3:!G 

- = - b<Mx') 

3z  3x'3z  3x'  3z2  3z:l 


We  substitute  again  for  x'  and  z  the  values  x  and  y  respectively  and 


write 

3G 

3G 

- -  u  and 

- —  — V 

3y 

dx 
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Hence,  we  obtain 


3u  3u  32u 

u - 1-  v - =  <f>  (x)  +  v - 

3x  3y  3y2 

and 

3u  3v 

- + - =  0 

3x  3y 

The  function  <f>  (x)  is  determined  from  the  initial  conditions, 
u  =  0  and  v  =  0  for  y  =  0;  therefore, 


or,  when  p  represents  the  pressure  distribution  along  the  wall,  then 
we  have  from  the  equation  of  motion  in  the  x  direction 

32u  \  1  3p 

-  )  - - =  —  <*>(*) 

3y2  /  y  =  0  p  3x 

If  this  value  is  substituted  for  t£(x),  we  obtain  the  differential 
equation  of  the  boundary  layer  which  was  given  by  Prandtl  in  1904. 

3u  3u  1  3p  32u 

u - h  v - = - 1-  v -  (6) 

3x  3y  p  3x  3y2 

This  equation  is  distinguished  from  the  first  equation  in  (1)  by 
the  omission  of  the  term  v  32u/3x2.  One  can,  of  course,  immediately 
understand  the  underlying  assumption  that  the  ratio  of  u  to  x  must 
be  vanishingly  small  compared  to  division  by  y;  in  addition  it  will 
be  found  that,  in  (1),  the  pressure  gradient  3p/3x  is  now  a  function 
of  x  alone.  It  can  be  proved  that  this  equation  is  not  changetl  when 
a  coordinate  system  is  introduced  at  a  point  on  a  smooth  curved 
wall  such  that  the  coordinate  distances  are  measured  normal  to  the 
wall  and  along  the  wall;  of  course,  it  is  assumed  that  the  curvature 
is  not  too  large.8  We  shall  always  use  this  coordinate  system  in  the 
future,  letting  the  origin  coincide  with  the  stagnation  point. 

The  thickness  of  the  boundary  layer,  8,  follows  from  the  con¬ 
sideration  that  it  is  proportional  to  a  specified  value  of  i jj  or  that  /  is 

*Sce  Hiemenz,  Dissertation,  p.  3. 
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of  the  order  of  magnitude  of  unity;  hence, 

8  — s  aR-1/2 

In  addition  we  note  that  the  Navier-Stok.es  equations  are  elliptic 
while  the  Prandtl  boundary  layer  equation  is  parabolic.  From  the 
physical  point  of  view,  the  Prandtl  equation  is  obtained  by  neglecting 
the  small  effects  within  the  boundary  layer  due  to  a  change  in  volume 
and  due  to  the  frictional  retardation  in  the  x  direction.  It  is  clear 
that  such  a  profound  change  in  the  differential  equation  alters  the 
calculated  streamlines  from  the  actually  observed  ones,  and  it  is  an 
additional  problem  to  determine  to  what  extent  the  flows  which  have 
been  calculated  by  boundary  layer  theory  will  agree  with  real  flows. 

The  important  phenomenon  of  separation  of  the  flow  from  a  sur¬ 
face  which  often  happens  when  fluid  flows  around  a  body  can  be 
explained  with  the  help  of  boundary  layer  theory.  When  the  thick¬ 
ness  of  the  boundary  layer  is  small  compared  to  the  size  of  the  body, 
the  pressure  gradient  across  a  cross-section  of  the  boundary  layer  is 
shown  to  be  nearly  constant.  The  velocity  distribution  in  the  bound¬ 
ary  layer  decreases  from  its  value  in  the  potential  flow  to  zero  at  the 
wall.  If  there  is  a  pressure  rise  along  the  body,  then  a  fluid  element 
within  the  boundary  layer  near  the  wail  will  come  to  rest  sooner  than 
one  located  in  the  potential  flow.  The  conditions  occurring  in  the 
neighborhood  of  the  separation  point  will  be  qualitatively  described 
by  Figure  1.  At  the  separation  point  itself,  one  has  the  following 
condition 

3u  32«|» 

_ — - =  0  at  y  =  0 

3y  3y2 

From  this,  the  position  of  the  separation  point  can  be  calculated. 
The  boundary  layer  separates  from  the  wall  at  this  point  and  moves 
at  a  definite  small  angle  into  the  main  stream  above. 

2.  Derivation  of  the  Karman  Integral  Condition. 

For  the  plane,  steady  state  problem  the  differential  equation  of 
the  boundary  layer  is 

3u  3u  32u  1  3p 

3x  3y  3y2  p  3x 

In  addition,  the  continuity  condition  is 
3u  3v 

- + - —  o 

3x  3y 
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We  define  a  function  q  (x,y)  such  that 
u  =  U  -  q  (x,y) 


(see  Fig.  2)  where  U  =  f,  is  the  velocity  outside  the  boundary  layer 
which  is  given  by  the  potential  flow.  This  potential  flow  is  derived 
from  the  combination  of  the  Eulerian  equation  and  the  experiment¬ 
ally  determined  pressure  gradient.  It  is 

1  dp  dU 

- =  U - =  UU' 

p  dx  dx 

Upon  introduction  of  this  function,  Equation  (6)  becomes 

/  dU  3q  \  3q  32q  til) 

(U-q) - -V - =  -v - +U - 

\  dx  3x  /  3y  3y2  dx 

or 

dU  3q  3q  3q  32q 

— q - U - (-  q - v - =  —  v - 

dx  3x  3x  3y  3y2 

We  integrate  with  respect  to  y  and  obtain 


The  upper  limit  may  be  written  as  x  provided  the  function  q 
goes  asymptotically  to  U.  We  can  assume  that,  at  some  distance  from 
the  wall,  the  thickness  of  the  boundary  layer,  q  does  not  vary  ap¬ 
preciably  from  zero.  It  is  enough  then  for  the  integral  to  extend 
from  zero  to  8. 

After  using  the  continuity  equation  and  after  integrating  by  parts, 
we  can  write 


x 

J3q  x 

v - dy  =  [vq]  - 

3y  0 

0 


S 


3v 

q  —  <*y 

3y 


=  [vq]  + 
0 


dU 

dx 


d 

dx 


dy 


0 


8 


00 

However,  [vq] 

0 

have  left 


vanishes  at  both  the  upper  and  lower  limits,  and  we 


v 


00 


dy 


(7) 


The  boundary  conditions  which  the  function  q  (x,y)  must  satisfy 
are 

1.  u  :=  0  i.e.,  q  =  II  at  y  =  0 

2.  u  =  U  i.e.,  q  =  0  at  y  =  oe 

3.  also,  the  boundary  layer  equation  must  be  applicable  at  the 
surface  of  the  body.  Since  both  u  and  v  are  zero  at  y  =  0,  it  follows 
that 


r  3*(U-q)  -] 

L  ay2 


y=0 


dU 

U - 

dx 


Every  solution  of  the  boundary  layer  equation  must  obviously 
fulfill  the  integral  condition  (7) .  If  we  now  make  a  simple  approxi¬ 
mation  for  the  dependence  of  q  on  y,  e.g.,  a  power  series,  the  coeffi¬ 
cients  will  be  only  functions  of  x.  An  approximation  to  the  solution 
is  obtained  if  one  determines  all  but  one  of  the  coefficients 
according  to  the  boundary  conditions.  The  integral  condition  then 
gives  an  ordinary  differential  equation  for  the  determination  of 
the  remaining  coefficient  and  thereby  fixes  the  dependence  of  the 
velocity  profiles  on  x. 


Although  the  original  equation  for  the  boundary  layer  is  a  partial 
differential  equation,  we  now  only  have  to  solve  an  ordinary  dif¬ 
ferential  equation  because  of  the  integral  condition.  It  is  not  difficult 
to  refine  the  approximation  further  by  satisfying  the -integral  con¬ 
dition  separately  for  several  subregions  (e.g.,  from  0  to  8/2  and  from 
8/2  to  8),  and  therefore,  instead  of  determining  one  remaining  co¬ 
efficient,  by  determining  two.  Also,  there  is  the  possibility  of  obtaining 
further  integral  conditions  in  the  following  way:  one  could  multiply 

both  sides  of  Eq.  (7)  by  y,  y2,  y:t . y",  one  after  the  other, 

and  could  integrate  from  zero  to  infinity,  and  then  perhaps  the 
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function  could  be  determined  in  analogy  with  the  Stieltjes  method 
(i.e.,  by  its  moments).* 

However,  the  following  calculations  show  that  the  fulfillment  of 
the  simple  integral  condition  preserves  the  physical  meaning  of 
boundary  layer  theory  and  determines  the  velocity  distribution  with 
sufficient  accuracy  for  practical  purposes.  The  particular  example  of 
the  flow  around  a  cylinder  shows  that  our  method  is  a  far  better 
method  of  calculation  than  a  previous  method  due  to  Hiemenz 
(Taylor  expansion  in  x  and  the  solution  of  a  series  of  ordinary 
differential  equations  for  the  determination  of  the  coefficients  of  the 
series)  since  after  the  solution  of  a  single  differential  equation  one 
obtains  values  which  are  comparable  in  accuracy  to  those  obtained  by 
the  solution  of  four  or  five  differential  equations  when  the  other 
method  is  used. 


It  is  shown  in  the  earlier  work  by  von  Karman”  that  the  integral 
condition  may  be  derived  on  the  basis  of  physical  considerations. 

3.  The  Floxv  Along  a  Flat  Plate. 


First,  we  will  use  our  approximate  method  for  the  integration 
of  the  differential  equation  of  the  boundary  layer  for  the  case  of  a 
flat  plate  when  the  plate  is  parallel  to  the  streamlines  of  the  uniform 
stream.  The  velocity,  U,  outside  the  boundary  layer  is  thus  a  con¬ 
stant,  and  the  differential  equations  (6)  read 

3u  3u  32u  3u  3v 

u - j-  v - =  v -  and - 1 - =  0 

3x  3y  3y'J  3x  3v 

The  integral  condition  (7)  simplifies  to 


3q 


8 


3y  _i() 


•Editor's  Note:  I’ohlhausen’s  suggestion  has  been  partially  exploited  by  l..  G. 
Loitsianskii,  XACA  Technical  Memorandum,  No.  1293,  1951,  who  used  the  zero 
(the  Ktirnuin  momentum  equation),  the  first,  and  the  second  order  moments  integral 
equations.  The  complete  Stieltjes  method  remains  to  lie  done. 


’'This  volume,  page  233  to  252.  An  English  translation  is  available  as  XACA 
Technical  Memorandum,  No.  1092.  1946. 
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The  boundary  conditions  are  u  =  0  at  y  =  0  and  u  =  U  =  const, 
at  y  =  oo  for  the  differential  equation  and  q  =  U  at  y  =  0  and 
q  =  0  at  y  =  oo  for  the  integral  condition.  In  addition 

(-)  =* 

\  3y2  /  \  =  o 

As  Blasius  showed  in  his  dissertation,  the  example  of  the  Hat 
plate  permits  the  reduction  of  the  partial  differential  equation  with 
the  help  of  a  similarity  transformation  to  an  ordinary  differential 
equation  of  the  third  order.  Blasius  integrated  the  continuity 
equation  by  use  of  the  stream  function,  tj/,  where 

d\ ft  3i ft 

u  = -  and  v  =z - 

dy  3x 

and  introduced  two  new  variables  f  and  £  in  such  a  way  that 


f 


Thence, 


1  /  pU  \% 

=  —  (  —  )  y  aml 

2  \  p  x  / 

p  \V4 

pxU  / 

1 

/  pU  ' 

u  =  %U£'  and  v  =  —  | 

17. 

1  (ft*  -  0 

where  the  primes  indicate  differentiation  with  respect  to  £.  If  these 
values  are  substituted  into  the  boundary  layer  equations,  we  obtain 


£"'  +  {£"  =  0 

The  numerical  integration  of  this  equation  has  been  worked  out 
by  C.  Topfer  by  use  of  the  Kutta  method."’ 

For  the  shear  stress  on  the  plate,  Blasius  obtained 


We  now  proceed  to  give  also  a  solution  for  the  Hat  plate  using 
the  approximate  method  based  on  the  integral  condition.  We  ex- 


niDescribed  by  C.  Runge,  /.eitscliroft  fur  Math,  unit  Phys.,  1912,  p.  597.  The 
development  of  the  function  £  and  its  derivative  is  also  given  by  K.  I’ohlhausen, 
I,  1921,  p.  119.  Fig.  1. 
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pand  q  in  a  power  series  in  y. 

u  =  U  —  q  =  a  (x)y  -f  b  (x)y!  -f  c  (x)yH  -j-  .  .  . 

1 .  The  series  is  terminated  after  the  first  term,  thus 
q  =  U  -  a  (x)y 

For  the  upper  limit  of  integration  we  insert  8.  We  also  replace  the 
boundary  layer  profile  by  a  curve.  The  boundary  condition  is 
u  =  U  for  y  =  8;  thence,  U  =  a8 

We  have  to  make  the  following  calculations: 


Therefore,  the  integral  condition  may  be  written* 
U28'  U28'  U  6V 

2  +  3  ~~  V  8  ~  U 


Consequently, 


For  the  shear  stress  we  then  obtain 


2.  Next,  we  obtain  a  better  value  by  approximating  the  velocity 
distribution  with  a  parabola  and  by  requiring  that  3q  3y  equals  zero 
at  y  —  8.  Of  course,  we  must  relax  the  condition  at  the  wall  that 
3 -’u  3y-’  =  0.  Thus  q  =  U  —  u  =  U  —  ay  —  by-’;  therefore,  a  and  b 


r  “18 

•Keillor's  Note:  The  value  of  _  must  he  taken  to  lie  — 11/8  because 

3y  j  () 

of  the  discontinuity  in  3<|/3v  at  \  =  g. 
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are  equal  to  2 U/8  and  —  U/82  respectively.  We  again  form  the 
necessary  integrals 


/US  n 

qdy  = -  anil  J  q2dy  - 


U28 


0 

and  it  follows  from  these  that 
S 


d 

dx 


f  -i-'i  = 


US' 


d 

dx 


r  U28' 

I  qady  = - 

J  5 


and 


3q 

3y 


0 

2U 


With  the  help  of  the  integral  condition  one  obtains 


U28' 


and  for  the  shear  stress 


1 

5 


2 pi!  /  SOpx  \ 

= -  and  8  =  (  -  ) 

S  \  U  ) 

b 

/  mPu  ’  V* 

0.365  (  -  J 


'I  he  first  approximation  yields  a  value  about  13  per  cent  too  small 
while  the  second  approximation  gives  a  value  about  6.5  per  cent  too 
large. 

3.  To  satisfy  both  conditions  (no  discontinuity  in  the  velocity  dis¬ 
tribution  and  tbe  boundary  condition  for  32q/3y2  at  the  wall),  at  least 
three  terms  must  be  retained  in  u.  We  let 

u  c=  ay  -)-  by-’  -j-  cy:‘ 

The  functions  a,  b  and  c  are  determined  from  the  boundary 
conditions: 

1.  u  =  U  for  y  =  8  gives  li  =  aS  -f-  b82  -f-  f8:i 

3u 

2.  - =  0  for  y  =  8  gives  0  =  a  -f-  2b8  -|-  3c82 


3y 


3. 


/  32u  \  UU' 

I  -  ]  = - =  0  and,  since  U'  —  0,  b  —  0. 

\  3y2  /  y  =  0  v 
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n 


We  obtain  then  the  equations  for  the  unknown  functions,  a  and 
c,  in  the  form: 

U  =  a8  +  c8s  and  0  =  a  -f-  3c82 
From  these  it  follows  that 


SU  U 

28  ’  283 


and  q  =  U 


3y 


The  integrals  become 

i  3U8  f 

J  qdy  =  — 'J  q2dy  = 

0  0 

and 


33U28 

140 


aq 


8  3U 


3y  [o  28 
The  resulting  differential  equation  for  8  is 

r  3  33  "1  $v 


us- 


140 


Therefore, 


28 


/  280j/x  /  vx  \,/4 

=  (-ss-)  =464  (-) 


y* 


28  28* 


The  friction  factor  for  this  approximation  is  0.323,  i.e.,  a  value 
which  is  about  3  per  cent  too  small. 

4.  Since  the  function  q  goes  to  zero  like  exp  (— tcy2)  according  to  the 
exact  boundary  layer  theory,  one  can  expect  that  a  still  better  approxi¬ 
mation  could  be  given  if  a  smooth  connection  to  u  =:  U  is  made  at 
y  -  8,  As  a  fourth  and  last  approximation  we  take  four  unknown 
functions,  a,  b,  c  and  d,  so  that 

u  —  ay  -f-  by2  -|-  cy''  -j-  dy< 
and  impose  the  following  additional  conditions: 

32u  UU' 

1.  for  y  =  0  :  - = - =  0 

3y!i  v 

2.  for  y  =  8  :  u  =  U 
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3.  for  y  =  S  :  3u/3y  =  0 
and  4.  for  y  =  8  :  32u/3y2  =  0 

From  these  we  see  that  the  equations  for  the  determination  of  the 
functions  are 

a  _|_  2b8  +  3c82  +  4d8<  =  0 
a8  4-  b82  +  c8:1  +  d84  =  U 
b  =  0 

and 

2b  -{-  6cS  +  12  d82  =  0 


The  solution  is 

a  =  2U/8,  b  =  0,  c  =  —  2U/83  and  d  =  U/S4 


so  that 


The  integral  condition  yields 

630y 

8'  = - 

37US 

and  from  this 


In  this  case  the  friction  factor  is  0.343.  This  value  is  about  3 
per  cent  too  large  so  that  this  approximation  has  about  the  same 
accuracy  as  the  preceding  one.  The  additional  calculations  required 
by  the  fourth  approximation  have  been  retained  since  the  smooth 
connection  of  the  velocity  distribution  to  the  external  How  at  y  =  8 


15 


lends  confidence  and  since  the  velocity  distribution  will  retain  its 
physically  sensible  form,  i.e.,  values  larger  than  U  will  be  avoided. 

The  values  of  the  friction  factors  which  have  been  given  by  the 
various  approximations  are  given  in  Figure  3.  In  order  to  obtain  a 
picture  of  the  development  of  the  boundary  layer  profiles,  the  cor¬ 
responding  profiles  are  determined  for  U  =  1,  x  =  1  and  v  =  1. 
(Compare  Fig.  4.) 


Approximation 


Exact 


y 

I 

II 

III 

IV 

Solution 

0 

0 

0 

0 

0 

0 

l 

0.289 

0.332 

0.318 

0.334 

0.329 

2 

0.578 

0.597 

0.606 

0.619 

0.629 

3 

0.867 

0.796 

0.834 

0.828 

0.846 

4 

1.000 

0.927 

0.972 

0.948 

0.955 

5 

1.000 

0.991 

0.990 

0.995 

0.990 

Finally,  let  the  "displacement  thickness” 

be  8*.  i.e., 

the  quantity 

indicating  the  outward  displacement  of  the  streamlines  of  the  po- 

tential  flow.  It  is  defined  by 


8 

1 

8*  = - 

U 

0 

For  the  various  approximations,  one  obtains 
Approximation 

1  II  III 

1.78  1.83  1.74 


IV  Exact 
Solution 


1.75  1.72 


4.  The  Flow  Around  «  Body  with  a  Given  Pressure  Distribution, 

The  original  idea  of  the  Prandtl  boundary  layer  theory  was 
that  Eq.  (6) ,  together  with  the  continuity  equation— in  combination 
with  the  theory  of  potential  How  of  an  ideal  Huid— sufficed  to  de¬ 
termine  the  velocity  distribution  in  the  boundary  layer.  Thus  the 
third  dependent  variable,  p,  which  appears  in  both  equations,  must 
be  eliminated  in  such  a  way  that  one  fixes  the  value  of  the  pressure 
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gradient,  dp/dx,  by  the  potential  solution  along  the  boundary  line; 
moreover,  the  value  IJ,  in  the  boundary  condition,  must  also  be  found 
from  the  potential  solution.  It  was  in  this  way  that  our  example  was 
treated  in  Section  3  with  dp/dx  =  0  and  U  =  const.  However,  the 
previously  mentioned  investigations  of  Blasius,  Boltze,  and  Hiemenz 
have  shown  that  no  satisfactory  agreement  with  observations  can  be 
found.  Hence,  following  a  suggestion  by  Prandtl,  Hiemenz  has  de¬ 
veloped  another  method:  he  experimentally  determined  the  pressure 
distribution  along  the  boundary  and  introduced  this  and  the  calcu¬ 
lated  distribution  of  U  into  the  differential  equation  and  the  bound¬ 
ary  conditions.  We  will  treat  the  case  of  a  pressure  distribution  given 
by  experiment.  In  this  general  case,  we  must  take  into  account  that 
separation  occurred  as  it  is  described  in  Figure  1.  Thus  it  is  useless 
to  represent  the  function  u,  which  must  describe  the  velocity  in  the 
boundary  layer,  by  a  polynomial  of  the  first  or  second  order  in  y 
since  the  first  algebraic  function  that  can  show  an  inflection  point  is 
one  of  the  third  order.  To  obtain  a  better  approximation,  we  will 
also  include  the  term  y4  and  will  write 

u  =  ay  -f  by2  -j—  cy*  +  dy4 

where  u  =  U  —  q  as  before.  The  thickness  of  the  boundary  layer  is 
again  8.  The  boundary  conditions  are 


for  y  =  6 

3u/3y  = 

for  y  =  8 

u  =  U 

92u 

for  y  =  0 

3y2 

for  y  =  8 

92u/3y2  : 

UU' 


where  the  fourth  condition  is  chosen  arbitrarily  as  we  did  in  the 
case  of  the  Hat  plate. 


The  four  unknown  functions,  a,  b,  c,  and  d,  are  now  determined 
with  the  aid  of  these  conditions.  From  condition  3.,  it  follows  im¬ 
mediately  that  the  function  b  can  be  determined  only  by  the  given 
pressure  distribution  and  is 

UU' 

~~  2v 
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The  remaining  conditions  furnish  equations  for  a,  c,  and  d, 
and  they  are 

a  +  2bS  +  3c82  +  4d83  =  0 
afi  -f  b82  +  c83  +  d84  =  U 

and 

2b  +  6c8  +  12d8'-'  =  0 

From  these,  the  pertinent  values  are  found  to  be 

U  UU'  UA 

a  - -  (12  +  A)  b= - = - 

68  2v  28- 

-U  (4  -  A)  U 

c  =  -  and  d  = -  (6  —  A) 

28‘  684 

where  A  is  the  dimensionless  number  U‘8'2/v. 

As  a  help  in  the  formulation  of  the  integral  condition,  we  now 
calculate 


0 


and  form 


0 


ant' 


3  8’ 

-f - U'8 - [  (U')2  +  UL’"] 

1(1  120«/ 


d  r  U-'8‘  p  11  U'8-  5  (U')284  “T  23 

-  1  q*dy  = -  23 - 1 - 7 - II  U'8 

dx  J  126  >_  4  v  12  v1  \  63 

0 

II  U8*  UU'8- 

- [2  (U')2  4-  UU"]  -f - [  (U') 2  +UU"] 

15 12  v  4536p2 
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•  ■•ft.**-  •  . .... 


(18 

US - 

dx 


If  we  put  these  expressions  into  the  integral  condition  (7)  and 
if  we  simplify  the  resulting  expression,  we  obtain  the  following 
differential  equation: 

llfi  84  U'8'1 

_2„  + _ l!'8- - [79  (U')2  +  HUH"] - [  (U')s  +  UU"] 

315  7560|/  4536i»2 

=  37  ui  5  (Uy84 

- 1 - 1 - 

315  315y  9072  v- 


or,  when  we  again  introduce  A  and  replace  8*/v  by  z. 


d/ 


dx 


LIU'1 


0.8  \  -9072  4-  I670.4A  -  47.4  +  4.8  • 


A" 


UU" 


1  4- 


(ur 


U  (-213.12+  5.76A+A-) 


The  solution  of  this  nonlinear  ordinary  differential  equation  ol 
the  first  order  gives  us  the  variation  of  the  thickness  ol  the  boundaiy 
layer  as  a  function  of  the  curvature  of  the  body. 


We  first  note  that  the  boundary  layer  has  a  definite  finite  thick¬ 
ness  at  the  stagnation  point  of  the  body;  thus,  the  differential  equation 
has  a  singular  point  at  x  =  0  and  8  =  8„  where  8,,  is  a  teal  positive 
root  of  the  cubic  equation 


r  uu"-j 

UU''- 

9072  +  1670.4 A  - 

47.4  +  4.8  1  A-  — 

(U?-l 

_l  +  ^ 

A‘  =  0 


To  determine  the  thickness  of  the  boundary  layer,  we  have  to  select 
the  integral  curve  which  passes  through  the  singular  jx)int.  The 
tangent  of  the  boundary  layer  profile  must  be  perpendicular  to  the 
wall  at  the  point  of  separation;  hence,  it  follows  that 

U 

a  = -  (A  +  12)  =  0  or  A  = -12 

68 

We  will  now  employ  the  universally  valid  differential  equation, 
which  we  previously  derived,  to  the  example  of  the  flow  around  a 
circular  cylinder  which  is  placed  in  a  uniform  stream.  The  pressure 
distribution  for  this  example  was  experimentally  determined  by 
Hienten/.  By  measuring  the  pressure  distribution,  his  experiment, 
which  was  performed  in  water  with  a  cylinder  of  97.5  mm  in  diameter 
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and  a  flow  velocity  in  the  undisturbed  fluid  of  19  cm/sec.,  gave 

U  =  7.15  lx  —  0.04497x‘  —  0.0003300x5 

where  the  ire  length,  x,  was  measured  in  cm.  We  form  U'  and  U" 
(see  Fig.  5)  and  assume  the  density  of  the  water  to  be  1  and  the 
viscosity  to  be  0.01. 

First  we  calculate  the  thickness  of  the  boundary  layer  at  the 
stagnation  point  as  the  root  of  the  cubic  equation;  this  calculation 
gives 

W 

A„  = - =  7.052  or  S„  =  0.09931  cm. 

v 

The  solution  of  the  differential  equation  can  be  obtained  either 
by  graphical  means  or  by  the  following  method:  substitute  a  series 
expansion  for  the  solution,  which  is  horizontal  at  the  singular  point, 
x  =  0  and  8-  v  =  h»/v,  and  then  extrapolate  this  curve  using  some 
method  of  calculation.  We  elect  to  solve  the  differential  equation 
graphically  by  use  of  the  Method  of  the  Isoclines.  Accordingly,  let 

dz  P  (x.z) 

dx  *  Q  (x.z) 

and  calculate  the  values  of  x  as  a  function  of  x  for  the  values 
z  =  8-  v  —  0,  1,2,  3,  4,  and  5.  These  values  are  obviously  infinite 
along  the  abscissae  for  which  the  denominator  vanishes.  This  is  the 
case  lor  A,  =  12  and  =  — 17.76.  Using  the  curves,  /  —  const.,  we 
now  construct  the  direction  field  for  the  differential  equation.  We 
draw  in  the  coordinate  system  x,  8-  v,  the  curves  \  =  const.,  and 
in  this  way  we  obtain  Figure  (i.  A  more  detailed  discussion  of  the 
differential  equation  reveals  that  the  singular  point,  x  =  0  and 
8-  v  —  8,,-V,  is  a  saddle  point."  Thus  only  two  isoclines— the  one 
with  a  horizontal  tangent  is  sought— pass  through  this  point.  If  we 
take  an  arbitrary  point  in  the  plane  :  id  attempt  to  construct  from 
this  point  the  integral  c  urves  by  use  of  the  isoclines,  then  the  t  urves 
return  to  either  plus  or  minus  infinity  because  of  the  above-mentioned 
singular  point.  In  this  way,  we  obtain  a  rigorous  criterion  for  the 
desired  integral  curve  p-  sing  through  the  singular  point.  The  arc 
length  to  the  separation  point  we  find  as  the  abscissa  of  the  inter- 

"See.  e.g.,  Horn,  ('•eu'iilinliche  /W/ZornfiVi/g/Wr/iiiNgcn  brlirliigfr  Oritnuiii’, 
l.cip/ig.  1903.  [>.  333. 
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section  ot  the  integral  curve  and  the  curve  X  =  —12.  This  gives 
x„  =  6.94  cm. 

Hiemenz  in  his  dissertation  has  made  an  attempt  to  solve  the 
differential  equation  of  the  boundary  layer  for  the  case  of  the  cylinder 
in  an  entirely  different  way.  With  the  aid  of  the  stream  function,  t \t, 
one  can  write 

j~  3 ip  3 ’■‘ip  difj  d-iji  ~  dp  d  if) 

p - - - K  - — 

3y  3x3y  3x  3ya _  tlx  3y:l 

He  expanded  the  pressure  gradient  dp  tlx  into  a  power  series  in  x. 
and  let  the  stream  function  be 

if/  =  i/t,x  +  t/*:,x:'  +  t//,x:' 

thus  describing  the  solution  in  the  form  of  a  Taylor  series  about  the 
origin.  The  partial  differential  equation  can  then  be  reduced  to  an 
infinite  scries  of  ordinary  differential  equations.  From  this  system, 
he  solved  the  first  three  without  testing  the  convergence  and  main¬ 
tained  that  tfj  is  determined  sufficiently  close  and  that  tj»7,  t/»„,  etc.,  do 
not  exert  a  substantial  influence.  He  obtained  a  value  ol  6.977  lor 
the  separation  point  which  is  approximately  the  value  derived  above 
and  which  agrees  very  well  with  experiment. 

As  a  result  of  graphical  manipulation  the  approximate  method 
has  a  great  advantage  in  giving  better  insight  into  the  development 
of  the  boundary  layer.  We  recognize  immediately  from  Figure  8, 
which  is  explained  below,  that  it  is  impossible  to  describe  correctly 
the  development  of  the  integral  curve  from  the  stagnation  point  to 
separation  point  with  three  terms  of  a  Taylor  series.  Also,  the  geo¬ 
metrical  location  of  the  separation  point — the  curve  X  =  —12 — re¬ 
sults  in  an  isocline  with  a  very  large  slope  so  that  considerable  change 
in  the  structure  of  the  integral  curve  brings  a  relatively  small  dis¬ 
placement  in  the  position  of  the  separation  point.  Hiemenz  thought 
that  this  loss  of  sensitivity  in  the  positi«:n  of  the  separation  point  was 
a  test  of  boundary  layer  theory. 

However,  the  difference  between  our  approximate  solution  and  the 
one  of  Hiemenz  was  very  considerable  when  one  calculated  the  pro¬ 
file  of  the  boundary  layer  for  some  value  of  x.  The  solution  ol  Hie¬ 
menz  altered  the  thickness  of  the  boundary  layer  only  very  little;  it 
amounted  to  about  1  mm. 

The  boundary  layer  of  the  approximation  agreed  at  the  stagna¬ 
tion  point  with  the  one  of  Hiemenz.  On  the  other  hand,  at  the 


21 


I 


separation  point  8*/v  —  384,  the  thickness  8  —  1.96  mm.  thus 
doubled  in  size. 

In  Figure  7,  the  displacement,  8*.  is  given,  both  according  to 
Hiemenz  and  according  to  the  above  solution.  It  is  defined  as 


o 

c  q  c 

q  ,  / 

y  \ 

8*  =  |  - dy  =  8 

—  d  ( 

J  u  J 

u  V 

s  / 

0  0 

If  we  substitute  the  ialue  U  —  ay  — by-  —  cy*  —  dy4  for  q  then 
we  obtain,  after  the  evaluation  of  the  integral  with  the  help  of  the 
dimentionless  parameter  X  =  U '82/v, 


The  following  table  gives  the  calculated  values  as  a  function  of  the 
arc  length  of  the  cylinder. 

x-cni  j()  1  2  3  4  5  5.5  6  6.5 _ 6^5  6^ 

8^v  ~  0.986  0.990  1.02  1 .08  1.21  L44  T62  h92  2T8  2^94  3.84 

8  mm  0.993  0.995  1.01  1.04  1.10  1.20  1.27  1.39  1.57  1.71  1.96 

8*inm  0.239  0.24  1  0.246  0.258  0.279  0.320  0.354  0.414  0.520  0.614  0.784 

We  suspect  that  our  solution  offers  a  better  approximation  to  the 
actual  development  of  the  boundary  layer.  In  order  to  show  the 
difference  clearly,  the  development  of  the  velocity  in  the  boundary 
layer  was  calculated  in  Figure  8  for  y  =  0.374  mm  (in  the  case  of 
Hiemenz,  H  =  I).  In  addition,  the  individual  terms  of  Hiemenz' 
series  development  have  been  shown.  This  shows  that  Hiemenz’  so¬ 
lution  probably  converges  to  our  approximate  solution,  that  the 
number  of  calculated  terms  if/,,  is  quite  insufficient  and  that  the  least 
one  should  do  is  calculate  t jt:  anil  if/,,.  It  appears  therefore  that  in 
this  case  a  series  expansion  about  the  stagnation  point  is  unsuccess¬ 
ful  since  the  variation  in  the  velocity  profiles  first  became  noticeable 
at  large  values  of  x  where  the  Taylor  series  converges  slowly. 

5.  The  Flow  Between  Non-Parallel  Walls  (Diffuser). 

An  additional  example  of  the  integration  of  the  boundary  layer 
equation  is  the  How  in  a  diffuser. 

We  treat  the  simplified  case  of  two  straight  and  nonparallel  walls 
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BC  and  DE  (Fig.  9) .  We  extend  the  straight  segments  to  the  inter¬ 
section  at  0,  and  let  the  distance  OB  be  unity.  OC  is  the  x  coordi¬ 
nate  in  our  system.  At  the  point  0,  we  assume  now  a  source  whose 
strength  is  2ttA.  It  is  then  obvious  that  the  radial  velocity  distribu¬ 
tion  is 

U  =  A./x 

This  value  is  now  put  into  the  universally  valid  differential  equation 
(8).  We  obtain  the  value  2  for  UU"/  (1)')-;  thus  the  differential 
equation  becomes 

dz  0.8  [—9072  +  1670.4  A  -57X-  — 3X3] 
dx  (A/x)[— 213.12 +  5.76X  +  X*] 

where 

dz  x 

X  —  zU'  —  —  Az/x-  and  - = - (2X  -)-  X'x) 

tlx  A 

If  these  values  are  inserted  into  the  differential  equation,  we  get 
x  dX  — 0.2X3  —  17.04X-  +  455.04X  —  3628.8 

2  dx  -213.12 +  5.76X  +  X- 

and  by  separation  of  variables 

2  dx  -213.12  +  5.76X  +  X-’ 

- - - dX 

5  x  18144  -  2275.2X  +  85.2X- -f  X3 

Separating  the  rational  function  on  the  right  hand  side  into 
partial  fractions  and  integrating,  we  note  that  two  of  the  terms  are 
zero  and  that  X  is  given  by 

2  107.85  +  X  X-  —  22.65X  +  168.2 

- lnx  —  0.7581  In - -f-  0.1210  In - 

5  107.85  168.2 

f  2X  -  22.65 

—0.0664  -I  arc  tan - (-  arc  tan  1.790 

[  12.65 

If  we  take  the  condition  for  separation  to  be  X  =  —12,  then 
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we  have 


2 

_ lnx„  =  0.7581  In  0.8887  +  0.1210  In  3.47 1 

5 

—  0.0664  [arc  tan  1.790  —  arc  tan  3.688] 
or  x0  =  1.214  for  the  separation  point. 

We  have  thus  obtained  the  remarkable  result  that  even  for  linear 
diffusers  with  an  arbitrarily  small  opening  angle,  the  boundary  layer 
theory  always  gives  a  separation  point;  the  latter  corresponding  to  a 
cross-section  widening  of  1.214:1,  i.e.,  about  21-22  per  cent.  However, 
this  result  is  not  physically  realistic  since  the  actual  pressure  distri¬ 
bution  deviates  from  the  one  based  on  the  present  analysis  in  the 
neighborhood  of  the  separation  point:  in  particular,  the  primary 
flow  is  mostly  turbulent. 

One  must  give  an  entirely  different  treatment  for  the  convergent 
channel.  We  place  a  sink  at  the  point  0  whose  strength  will  be  de¬ 
noted  again  as  i.  For  this  case,  the  radial  velocity  distribution  is, 

V  =  —  A/x,  and  the  differential  equation  is  now 

x  dX  X1  +  I2.98X-  -  400.6X  +  1649 
T?"dT_  213.12  —  5.76X  —  Xs 

The  integration  proceeds  again  by  the  method  of  separation  of 
variables 

dx  213.12  -  5.76X  -  X- 

4.4 - - - dX 

x  X:'  +  I2.98X-  -  400.6X  +  1649 

and  yields 

,1  Inx  =  0.2173  In  (X  -  10.34)  -  0.91 15  In  (X  -  5.530)  -  0.33631n  (X  +  28.85)  +  C 
As  anticipated,  this  flow  gives  no  separation  point. 

The  rigorous  integration  of  the  Navier-Stokes  equations  has  been 
accomplished  by  others  for  the  previous  example  of  the  sink  with 
infinitely  long  walls,  thus  a  comparison  should  be  made  between 
their  solution  and  the  approximate  solution.* 


•Kilitor’s  Note:  The  suggested  comparison  was  made  try  K.  Millsaps  ami  k. 
I ‘oh  I  ha  usen.  lour.  Aero.  Sri.,  Vol.  20.  p.  191.  I9.V1 
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We  pul  the  stream  function  in  the  form 
ip  —  ip  (t)  where  t  =  y  x 


and  obtain 


dip  ip' 


dip  yi p‘ 


dx  x- 


du  yip"  ip'  3u  ip"  9-u 

dx  x-‘  x-  3y  x-  dy~ 

The  pressure  gradient,  dp  /dx,  is  calculated  as 

I  dp  A- 

- -  mi'  = - 


If  we  substitute  these  values  into  the  boundary  layer  equation,  we 
get  an  ordinary  differential  of  the  third  order 

—  WY  =  —A8  +  i ip'" 

Letting  ip'  equal  £, 

-£8  =  —A-  +  v{" 


-*Af  =  -A 8{  4- 4  i;  (O2  +  » 

The  integration  constant  B  is  determined  by  the  condition  that 
£  —  — \  for  £'  =  0;  it  then  follows  that 

2  v  £  2A‘ 

B  = - A:t  and -  (£')-  = - f-  A-£  -j - 

3  2  3  3 


(^r 


i'  =  (£  +  A)  (2A  —  £) 


(-9  ■  (4)' « 


+  A)  (2 A  -  £),/8 


We  now  make  the  substitution  q  =  2  (2A  —  £)’  - 
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or 


dq  = - 

(2A  -  £)•/» 

and  then  obtain 


ISO 

> 

1 

(  A  \‘i  y  - 

)  =  tanh 

a+  (  -  )  - 

\  3A  , 

! 

L  V  2  v  )  x  _j 

The  integration  constant  a  is  determined  from  the  condition, 
£  =  0  for  y  =  0, 

2  —  3tanh2a  =  0,  tanh  a  =  ^ ^  or  a  —  1.146 
so  that 


£  =  2A  —  3A  tanh2 

The  velocity  profiles  for  the  exact  solution,  for  A  =  1,  v  =  1. 
x  =  1,  and  for  the  approximate  solution  are  shown  in  Figure  10. 
The  root  which  is  used  for  the  approximate  solution  is  10.34. 


1.146  -f 


/  A  \  Mi  y 
\  2v  )  x 
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Figure  2.— The  Velocity  Deficiency  Function. 
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Figure  3.— Friction  Factors  for  a  Flat  Plate. 
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Figure  4.— Velocity  Profiles  at  a  Point  on  a  Flat  Plate. 


Figure  5.— Velocity  of  Flow  Outside  the  Boundary  Layer  on  a  Circular  Cylinder. 
Experimental  Data  by  K.  Hiemenz. 
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Figure  10.—' Velocity  Profiles  in  a  Converging  Channel  According  to  Boundary  Layer 

Theory. 
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